Nonlinear Transport through NS Junctions due to Imperfect Andreev Reflection 
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We investigate a normal metal - superconductor (point) contact in the limit where the number 
of conducting channels in the metallic wire is reduced to few channels. As the effective Fermi 
energy drops below the gap energy, a conducting band with a width twice the Fermi energy is 
formed. Depending on the mode of operation, the conduction band can be further squeezed, leading 
to various non-linear effects in the current-voltage characteristics such as current saturation, a re- 
shaped negative differential resistance, bistability, and hysteresis. 
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Coherent transport phenomena in micro-structured 
normal-superconductor (NS) systems have recently at- 
tracted a lot of interest M. The transport across a (NS) 
boundary is governed by the phenomenon of Andreev re- 
flection [p| : An electron incident from the normal metal 
on the NS junction with an excitation energy |e| below 
the superconducting gap A cannot penetrate into the 
bulk superconductor (we measure the excitation energy 
e of the electron with respect to the chemical potential fi 
in the superconductor). Nevertheless, subgap transport 
across the junction is possible via the process of Andreev 
reflection, where the electron incident on the boundary 
is accompanied by a (coherent reflected) hole, produc- 
ing effectively a state with two incoming electrons which 
convert into a Cooper pair upon entering the supercon- 
ductor. For an ideal NS boundary, such a process leads to 
a conductance G = Ae 2 /h per channel ||, twice higher 
than the maximally possible normal one. If the trans- 
parency T of the boundary is smaller than unity, the NS 
linear conductance decreases as T 2 at small T < 1 ||. 
New effects appear in the finite-bias or finite-temperature 
conductance when the transmission of electrons and holes 
differs significantly as a consequence of their different lon- 
gitudinal kinetic energies §J|-§]. In this letter we show 
how Andreev scattering, combined with specific condi- 
tions for the propagation of electrons and holes, leads to 
the formation of a subgap conduction band with a width 
which strongly depends on the bias voltage, leading to 
new transport characteristics exhibiting a negative dif- 
ferential resistance, bistable, and hysteretic effects. 

To fix ideas, consider a single channel normal metal 
point contact to a bulk superconductor, see Fig. 1(a). 
Here we have in mind geometries such as quantum point 
contacts realized in heterostructures j^] or via manipu- 
lations with a scanning tunneling microscope Given 
the chemical potential /i, we define the longitudinal chem- 
ical potential \i x = ji — £j_, where e± denotes the trans- 
verse energy of quantization in the normal channel. Here, 
we are interested in situations where the longitudinal ki- 



netic energy K e ^ at the Fermi surface is small, such that 
the condition /j, x < A is realized. In this case, electrons 
with excitation energies e = E — ji = K e — ji x > can 
propagate through the contact, whereas the correspond- 
ing hole state with the same excitation energy can prop- 
agate only if its kinetic energy Kh = /j. x — £ remains 
positive, i.e., e < fi x ( m this simplified discussion we as- 
sume that the conducting normal channel is long enough 
to generate transmission probabilities or 1 only) . 

An electron incident on the superconductor defines the 
Andreev state (in the normal single channel region) 
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with k± = y / 2m(/i 2; ± e)/h and v± = Tik±/m [the above 
states are normalized to carry unit particle flux, with a 
normalization (^ £ ,^ e /) — 2tt?iS(e — e'), implying that 
|r ee | 2 = 1 — |r e ft | 2 ] . Following the above discussion, the 
quenching of the hole state for energies e > \i x , combined 
with the restriction in the allowed energies for incident 
electron states e > —/J, x , leads to the formation of a con- 
ducting band of width 2fi x , see Fig. 1(b) and the inset of 
Fig. 2. Within this band, incident electrons are (nearly) 
perfectly reflected into holes, whereas electrons with en- 
ergies above this band (e > fi x ) are reflected as elec- 
trons and do not carry current (electrons with energies 
(e < —fi x ) do not enter the normal channel at all). 

In the simplest formulation of the problem we consider 
a single channel NS junction. The Andreev states are 
found by solving the Bogoliubov-de Gennes equations 
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with the gap function A(x) = A0(x), using the Ansatz 
(|l|) in the normal region and 
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in the superconducting region x > 0, with p 2 — q 2 = 
2mp x /h 2 , pq = mVA 2 -e 2 /h 2 , 7 = (e- WA 2 ~e 2 )/A, 
and the normalizing velocity v — hp/m. Solving for 
the transmission and reflection coefficients we obtain the 
spectral conductance Gns(|e| < A) rss (8e 2 /h)®(fi x — 
\ £ \) \fT^~~^ I A, valid in the limit fi x -C A (see in- 
set of Fig. 2; note that the Andreev approximation is 
not valid in the present case). Assuming a rigid band, 
implying that the applied bias drops at the boundary 
of the reservoir to the normal lead, the finite conduc- 
tion band produces a current-voltage (I-V) characteris- 
tics I — J c de G*Ns/e which saturates at a bias eV = fi x 
(we assume a negative bias V < 0, hence eV > 0). This 
then is the simplest example where the quenching of the 
back-propagating holes limits the width of the conducting 
band and entails a non-trivial saturation phenomenon in 
the transport characteristics of the NS junction. It is in 
contrast to the normal point contact, where an increasing 
bias eV opens up the channels, see Refs. |nj (note that a 
non-trivial structure in the I-V characteristics can occur 
in a normal point contact as well, though at much higher 
voltages eV > see PH). 

The above ID example cannot be trivially applied to 
a physical situation: strictly speaking, we neither have a 
superconducting instability nor a normal Fermi liquid in 
a single channel ID system. Let us then discuss the more 
realistic single channel NS point contact with a geometry 
as sketched in Fig. 1(a), see also 0: The ID normal wire 
is adiabatically connecting the normal reservoir with the 
bulk superconductor. The smooth form of the wire guar- 
antees the appropriate matching of the wavefunctions in 
the three device segments, reservoir, ID wire, and bulk 
superconductor. An imperfect matching leads to a nor- 
mal reflection of the incident particle at the NS boundary 
and thus reduces the subgap transport with its interest- 
ing new features (below we will discuss the consequences 
of normal reflection for the device characteristics in more 
detail). The main properties of this geometry are the 
following: i) the confining potential of the wire produces 
a narrow conduction band connecting the reservoir with 
the NS boundary, ii) the ID wire is long enough to guar- 
antee a sharp onset of the transmission (but short enough 
to let us ignore strong interaction effects due to the one- 
dimensionality), iii) the wide NS contact helps the proper 
matching of the wavefunctions, and iv) the large chemical 
potential in the bulk superconductor allows us to adopt 
the Andreev approximation. The functionality of this 
device resembles that of the idealized structure above: 
whereas in the ID wire the back-propagation of the hole 
was limited by the bottom of the conduction band, in the 
present situation holes reflected from the NS boundary 
with minimal kinetic energy \x x — e < have to tun- 
nel through the effective potential due to the transverse 
quantization in the wire, see Fig. 1(c). 

As we can apply the Andreev approximation for the 



geometry of Fig. 1(a), the determination of the conduc- 
tance is trivial, G NS (e) = (ie 2 /h)Q(p x — \e\), and the cor- 
responding I-V characteristics for the rigid band model 
follows from simple integration, see Fig. 2. 

So far, the determination of the current-voltage char- 
acteristics has been based on a rigid band model, where 
the voltage drop in the device occurs at the boundary 
to the normal reservoir. A more accurate calculation of 
the transport current I(V) involves a self-consistent de- 
termination of the charge p(x) and the electric potential 
e<p(x) in the wire, given an applied bias V, see, e.g., Ref. 
p2| . Here we refrain from such a calculation, but rather 
discuss two interesting limiting cases illustrating the po- 
tential features of such a device. 

The first case we wish to analyze is the gated wire, 
where a top gate placed over the wire modifies its po- 
tential eSip = U(V, V g ). In the simplest case, considered 
by Brown et al. O], the gate potential V g follows the 
applied bias, 5V g — V . Assuming further that the wire 
potential is slaved to the gate, the band bottom in the 
wire is lifted by U — eV, implying that backpropagat- 
ing holes and low energy incident electrons are cutoff at 
|e| = ±(/i x — eV) rather than ±/x x . Within the Andreev 
approximation the spectral conductance G NS is narrowed 
down to the interval |e| < p x — eV and takes the form 



G ns (s,U(V)) 
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- U(V)) - \e\}. 
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A simple integration produces the I-V characteristics 

I(V < A/e) = - f deG NS (e, U{V)) (5) 
e Jo 

exhibiting a negative differential resistance (NDR) re- 
gime within the bias interval p x /2 < eV < p x , see Fig. 2: 
Increasing the applied biased eV up to p x /2, additional 
current carrying states are occupied and the transport 
current / increases. Going beyond the value p x /2, the 
rising bottom of the band quenches the back propaga- 
tion of the holes and fewer states are available, until at 
eV = p x all current carrying states are blocked. 

The above NDR phenomenon is further accentuated 
in a device where the charge of the wire rather than its 
potential is fixed — this is the second limiting case we 
wish to study here. The contribution to the charge den- 
sity of an individual channel averaged over the wire cross 
section is given by 

p(x) = 2e]T{/ v ( £fc )Mz)| 2 + [1-Me k )]\v k (x)\ 2 }, (6) 

k 

with f v (e) the (bias dependent) distribution function for 
the Bogoliubov quasi-particles. We evaluate the density 
in the middle of the wire and allow for a non-zero poten- 
tial shift U. Using the normalization introduced in (|l|) 
we arrive at the form 



P= ^h l d£ 
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with v± — \j2{p x — U ± ej/m the velocities of the quasi- 
particles. For the case of perfect Andreev reflection the 
above expression simplifies to [we assume an open chan- 
nel configuration with < eV < p x — V '; the occupation 
numbers are determined by those of the metallic reser- 
voir, / v (e) = 0(— e + eV) at zero temperature] 
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Requiring that the charge difference Sp = p(eV, U) — 
p(0, 0) vanishes at any applied bias V leads to the condi- 
tion (k F ^ x = y/2mp x /Ti) 



ek F x 



Hx 



0, 



determining the potential shift U {V) in the wire. Solving 
for U we obtain the result 



U(V) 



{eVf 

^Px 



< eV < p x - U{V). 
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The negative shift in the wire potential seems quite puz- 
zling at first sight, but can be easily understood in terms 
of the reduced group velocity of the back-propagating 
holes. A second solution is found for a positive shift U > 
p x , where the Andreev scattering is quenched and all inci- 
dent electrons are reflected back as electrons (l^eel 2 = 1)- 
With Sp = (2efc F ^/7r)[^l - (U - eV)/p x - 1] = we 
find the shift 



U(V) = eV, p x - U(V) < < eV. 



(10) 



Finally, a regime with partially quenched Andreev scat- 
tering is found in the interval < p x — U < eV, 
with U(V) determined by the relation dp = (ek FX /7r) 
[2^\-{U ~eV)/p x - ^2{1 - U/p x ) - 2} = and the 
result 



U = 2eV - 5p x + A^/p x {2p x - eV). 



(11) 



The internal potential shift U versus applied bias V is 
shown in Fig. 3 (thick solid line). The three branches 
of U(V) exhibiting completely open, partially quenched, 
and entirely quenched hole propagation in the wire ar- 
range to define a typical bistable configuration of the 
wire within the bias interval p x < eV < 2p x . The lower 
branch with the open channel terminates at eV = 2p x 
and the system has to jump to the state where the back- 
propagation of holes is quenched. Physically, the jump 
between the two branches corresponds to a rearrange- 
ment of the potential drop in the device: At small bias 
(lower branch) the applied bias drops on the left side of 



the channel, towards the reservoir. At high bias (upper 
branch) the potential drops on the superconductor side, 
producing the gated situation described above. Translat- 
ing this behavior of the internal device bias U to the I-V 
characteristics, see Fig. 2, we find a jump from the 'open' 
current carrying state at low bias to the 'closed' gated 
state at high bias as the applied bias eV grows beyond 
the band width 2p x , thereby producing a characteristics 
with a N-shaped negative differential resistance. At volt- 
ages eV > A a finite conductance is restored. Note that 
the I-V characteristics is not symmetric: for a positive 
applied bias V > (eV < 0) the conduction band stays 
open up to the bias eV — —2p x , where the wire potential 
aligns with the potential in the reservoir, U — p x = eV, 
see Increasing the bias further, the current saturates 
(similarly to the rigid band case) as part of the incoming 
electrons are excluded from entering the wire. 

The above analysis for the case of ideal Andreev reflec- 
tion at the NS boundary can be easily generalized to take 
a finite normal reflectivity of the barrier into account. 
It is convenient to characterize the junction through its 
normal state properties: Given the reflection coefficient 
R for electrons entering the ID-wire, the parameter \r ee \ 2 
switches between the values \r ee \ 2 = AR/ (1 + R) 2 ('open' 
channel) and \r ee \ 2 — ('closed' channel) (see Refs. 
fn|0|; we assume that d e R(e) w and e <C A, allowing 
us to ignore the energy dependence in |r ee | 2 ). A finite 
value of R then leads to a smoothing of the potential- 
bias relation and the current-voltage characteristics, see 
Figs. 2 and 3. As R approaches unity the fixed-charge 
characteristics approaches the result of the gated wire. 

Above we have concentrated on the single channel 
limit, where the non-linear effects leading to the NDR 
phenomenon are most pronounced. Going over to a 
multi-channel geometry, our analysis can be carried over 
to the grazing incidence trajectories Q with the modifi- 
cation p x — > p Xin = p — s± >n , where e±^ n = h 2 K 2 /2m 
is the transverse energy of the n-th channel. The in- 
teresting structure obtained in the single channel case 
(see Fig. 2) then survives for the states with an effective 
chemical potential p x ,n < A. The maximal number of 
channels to be saturated or cut off in this fashion below 
a bias eV ~ A is of order Sn ~ (A//i)n, where n S> 1 
is the total number of channels. In a bulk system with a 
planar NS boundary these states correspond to grazing 
incidence trajectories with angles $ < y/A/p. 

In conclusion, the different propagation conditions for 
electrons and holes in NS junctions produce a rich variety 
of phenomena. Both, the zero- and finite-bias anomalies 
[ 1 1 in dirty NS junctions can be understood in such terms 
[S|. Here, we have shown how the coherent electron-hole 
transport in a one/few channel system may lead to strong 
non-linearities in the device characteristics, resulting in 
an N-shaped NDR I-V curve in its most extreme vari- 
ant. In conventional semiconductor devices this kind of 
instability leads to the formation of domain walls, e.g., 
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the Gunn effect. In the present case where the transport 
is non-local and coherent we can expect a device oper- 
ation more similar to that of a double-barrier resonant- 
tunneling structure | OL 
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FIG. 2. Current-voltage characteristics of the NS contact. 
Keeping the band rigid (dashed line), the current saturates 
when all electrons within the band 2fi x contribute. In the 
gated wire (long dashes), the back-propagation of holes is 
partially inhibited when the voltage increases beyond fi x /2 
and completely quenched beyond eV = fi x . In the wire with 
fixed charge (thick solid line) the current switches from the 
upper to the lower branch, producing a pronounced negative 
differential resistivity. A finite reflectivity reduces the insta- 
bility (fixed charge, thin solid lines, \r ee \ 2 — 0.25, 0.5, 0.8). 
Inset: Spectral conductance G N s versus energy e. For the 
ID NS wire the conductance is suppressed due to imperfect 
Andreev reflection (dotted line). For the adiabatically con- 
nected wire of Fig. 1(a) the Andreev approximation is appli- 
cable and Gns reaches is maximal value. The width of the 
conduction band depends on the wire potential U (solid line: 
gated/fixed-charge wire; dashed line: rigid band). 



FIG. 1. (a) Geometry of the NS point contact with the ID 
normal wire adiabatically connected to the normal reservoir 
on the left and the bulk superconductor on the right, see Ref. 
[8]. (b) Dispersion relation in the normal wire. Note the for- 
mation of a conducting band (cb) of width 2/j, x . (c) Energy 
diagram for the wire sketched in (a). A change in bias eV 
induces a shift U in the wire potential, which in turn may 
lead to a reflection of the back-propagating hole. The elec- 
tron then is normal-reflected from the NS boundary and does 
not contribute to the current. 



FIG. 3. Potential U within the normal wire versus ap- 
plied bias for the case of a fixed charge and different val- 
ues of the reflection coefficient r ee , r ee = (thick solid line), 
\r ee \ 2 = 0.25, 0.5, 0.8 (thin solid lines). The lower branch 
becomes unstable at high applied bias and the internal po- 
tential U jumps to the gated value (long dashes), leading to 
a NDR in the I-V characteristics of the contact. 



4 




(a) 



reservoir ID-wire 



super 
conductor 




(c) 



super 
conductor 



j'eV 




U 

ID-wire 




